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NOTATION 


a  Radius  of  cylinder  In  z-plane 

a^  Coefficient:  in  Laurent  series  (Bl) 

b  Distance  from  body  center  to  fin  tip 

bq  Coefficient  in  Laurent  Series  (B2) 

B  Function  in  Equation  (13) 

B^  Function  in  Equation  (14) 

c  Radius  of  circle  in  C-plane 

k  Coefficient  in  Equation  (4) 

Total  strength  of  the  point  vortices  behind  the  l^-fin 

l  Coefficients  in  Equations  (3)  and  (9) 

m  Total  number  of  point  vortices  in  the  field 

M  Number  of  point  vortices  representing  a  single  discontinuity  line 

n  Number  of  fins 

p,q  Coefficients  in  Equation  (7) 

r,<t>  Polar  coordinates  (r  also  coefficient  in  Equation  (5))  • 

Re  Reynolds  number  *  2bU/v 
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s  Coefficient  in  Equation  (6) 
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w  Complex  potential 

Wp  Contribution  to  w  From  parallel  flow 
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PREFACE 


This  report  Is  part  of  a  continuing  effort  at  the  Computation, 
Mathematics,  and  Logistics  Department,  with  support  from  IR-inhouse 
funds  and  the  6.1  NAVSEA  Mathematical  Science  Program,  to  study 
vortex  shedding  from  solid  bodies  in  a  fluid  flow  and  to  apply  the 
results  to  Navy  problems.  For  the  last  ten  years  the  major  objective  of 
this  effort  has  been  to  investigate  vortex  generation  and  shedding  in 
real  fluids  by  the  numerical  solution  of  the  Navier-Stokes  equations. 
These  successful  studies,  which  resulted  in  numerous  publications  in 
the  open  literature,  were  originally  restricted  to  moderate  Reynolds- 
number  flows  about  simply-shaped  bodies.  Today  two-dimensional 
flows  around  bodies  of  quite  arbitrary  shape  can  be  handled,  but  the 
solution  of  the  Navier-Stokes  equations  for  high  Reynolds  numbers 
still  cannot  be  obtained.  Instead,  ideal  fluid  flow  models  with  their 
well-known  shortcomings  must  be  used.  This  report  describes  one  of 
two  preliminary  studies  to  develop  a  computer  program  for  vortex 
shedding  past  arbitrarily  shaped  cylindrical  bodies  within  the  realm  of 
ideal-fluid  models.  This  report  deals  with  vortex  shedding  from 
finned  cylinders,  and  the  forthcoming  second  report  by  R.  Shonff  will 
address  vortex  shedding  from  arbitrarily  shaped  bodies  excluding  fins 
and  other  sharp  protuberances.  These  purely  two-dimensional  flows 
then  may  be  used  in  a  strip  theory  to  Include  at  least  some  aspects  of 
three-dimensional  flows.  The  ultimate  goal  will  be  a  computer  code 
for  vortex  shedding  from  three-dimensional  bodies. 
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ABSTRACT 


l 


A  computer  program  has  been  developed  to  simulate 
vortex  shedding .from  circular  cylinders  with  evenly  distri¬ 
buted  fins.  These  bodies  may  rotate  continuously  or  oscillate 
in  a  parallel  stream.  The  tyb-dimenslona]./flow  model  consists 
of  point  vortices  inserted  in  an  otherwise  potential-flow 
field.  For  the  roll-up  o,f  the  line  of/vortices  the  rediscre¬ 
tization  scheme  by  Fink/£nd  Soh  is  used.  Sample  results  are 
presented  for  vortex  shedding  from  /  flat  plate  at  angles  of 
attack  of  45®  and  90  ,  and  from  a/dircular  cylinder  with  two 
fins  at  an  angle  of  attack  of  45  .  . 

1  ~  . '  / 

ADMINISTRATIVE  INFORMATION  ( 


The  work  presented  in  this  report  was  supported  by  the  Independent  Research 
Program  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  under 
Work  Unit  1843-050,  and  the  6.1  NAVSEA  Mathematical  Sciences  Program  under  Work 
Unit  1808-010. 


1.  INTRODUCTION 

The  simulation  of  vortex  shedding  from  bodies  in  potential  flow  by  means  of 
point-vortex  models  has  attracted  the  attention  of  many  researchers  in  the  last 
decade  for  a  number  of  reasons.  Persisting  difficulties  in  solving  the  Navier- 
Stokes  equations  tor  large  Reynolds  numbers,  the  availability  of  large  computers, 
and  progress  in  the  study  of  rolled-up  discontinuity  sheetB  have  foBtered  the  use 
of  point-vortex  models.  The  extensive  literature  on  this  subject  includes  recent 
survey  papers  by  Fink  and  Soh,1*  Saffman  and  Baker,2  Clements  and  Maull,3  Koto,1'  and 
Leonard . b 

Although  the  neglect  of  viscosity  limits  the  usefulness  of  point-vortex 
methods,  in  many  cases  details  of  the  flow  field  can  be  obtained  and  a  fairly  good 
estimate  of  the  force  coefficients  can  be  made. 

This  report  presents  the  equations  of  motion  for  incompressible  fluid  flows 
past  abruptly  started  circular  cylinders  with  n  evenly  distributed  fins  of  equal 
length.  These  cylinders  may  rotate  continuously  or  they  may  oscillate.  Some  cases 
can  be  extended  to  elliptic  cylinders.  Point  vortices  are  introduced  into  the 
potential  flow  around  such  cylinders  to  simulate  the  development  and  shedding  of 


*A  complete  listing  of  references  is  given  on  page  39. 
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vortices  at  the  fins.  A  fixed  interval  of  time  elapses  between  successive  intro¬ 
ductions  of  point  vortices.  The  computer  program  is  checked  for  simple  cases  by 
comparing  its  results  with  solutions  from  the  literature.  Results  for  force 
coefficients  and  complicated  flow  problems  of  practical  interest  will  be  presented 
at  a  later  time  in  another  paper. 

The  formulas  derived  in  this  report  and  the  computer  program  described  here 
can  be  applied  in  missile  aerodynamics  and  ship  hydrodynamics.  In  particular,  cross 
flows  past  cruciform  fin  configurations  and  past  underwater  vehicles  with  sails, 
rudders,  stabilizers,  bilge  keels,  and  cables  can  be  determined. 


2 .  CONFORMAL  TRANSFORMATIONS 


For  the  computation  of  the  flow  field  the  method  of  mapping  the  physical  plane 
z  -  x  +  iy  onto  the  circle  plane  ?  ■  £  +  in  is  used.  Numerical  methods6*7*8  are 
available  which  map  an  arbitrarily  shaped  body  contour  onto  a  circle  by  means  of 
the  transformation  z  ■  f(?).  While  our  work  was  in  progress,  Mendenhall,  Spangler, 
and  Perkins8  published  a  paper  on  vortex  shedding  from  arbitrarily  shaped  bodies 
using  a  numerical  mapping  technique  for  the  Theodorsen  transformation.  V.A.  Golovkin 
and  M.A.  Golovkin10  worked  with  Fredholm  integral  equations  to  compute  the  roll-up 
of  point  vortices.  In  this  report  exact  conformal  transformations  f(c)  are  applied 
to  avoid  errors  due  to  the  approximation  of  the  body  contour.  Of  course,  these 
exact  transformations  are  restricted  to  certain  classes  of  bodies.  For  a  cylinder 
with  a  circular  cross-section  of  radius  a  and  with  n  evenly  distributed  fins  of 
length  b-a,  Miles11  has  given  f(r.)  in  the  implicit  form  (eee  Figure  1) 


+  <a2/z)(n/2)  -  c(n/2)  +  (c2/0(n/2) 


2c(n/2)  -  b(n/2)  +  (u2/b)(n/2) 


(1) 

(2) 


where  c  is  the  radluB  of  the  circle  in  the  c-plane.  Equation  (1)  crn  be  written 
in  the  explicit  form 


2  n  n  rff - n - 12 

z  -  2  n  U2  +  <c2/r.)2  +  /[  C2  +  <c2/c>2  ]2  -  4an}n  (3) 

In  general,  this  expression  is  not  single-valued  and  care  must  be  taken  in  working 
with  it.  From  Equation  (3)  it  can  be  shown  that  (dz/dc)  -  1. 
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For  fina  protruding  in  a  dihedral  configuration  at  an  angle  different  from 
90°,  an  approximate  conformal  mapping  is  given  for  n-2  in  implicit  form  with  an 
auxiliary  mapping  in  the  Zg-plane  (Figure  2) : 


z  +  —  ■  z,  -  - - (b  +  £-)  coe  0 

Z  L  *2  D 


Z2  + 


r2  c2  i 

-<+r+8+f 


(4) 


(5) 


with 


k  "  |(b‘b~)sin  0 
fl,  ■  ^(b  -  a)2  cos  0 


r  -  /|r t  +  4k2) 

...|(A3  +  B3),  c.i(A3-E3) 


E 


2  E 
k 


-2a  +  (b  +^-)sin  0  for  E2  <  0 


A2  -  X~  m  2a  +  (b+  §“)sin  0  for  A2  >  0 


A  ■  A  +  — — — t —  -  -- 
"3  2  a2-S./2  2 


E3  '  E2  +  ipi7I  -  f  «> 

Here  again,  (dz/dC)  ■  1.  The  derivation  of  this  conformal  transformation  is 
£  ■0° 

given  in  Appendix  A. 

In  certain  cases,  for  instance  for  two  fins,  the  circular  cylindrical  body  can 
be  replaced  by  an  elliptic  cylinder: 


where 


z  +  -£p+a)2  _  t  m  r  +  si 

1  4z1  v  4; 


c.|tTi  +  ^  +  l2  +  toil, 


tj  -  fttj  +  /  tj  -  pS  +  q2'  ) 


t2  -  i(t2  +  / 1\  -  P2  +  q2  )  (9) 

(see  figure  3),  Equations  (7)  through  (9)  are  Bimilar  to  those  given  by  Bryson.12 
However,  Bryson's  formulae  are  not  entirely  correct,  &a  shown  by  the  case  a»b, 
tj^-tjj-a,  s  arbitrary  (in  Bryson's  notation). 

3.  FLOW  FIELD 

3.1  COMPLEX  POTENTIAL  AND  COMPLEX  VELOCITY 

In  the  problem  under  consideration  the  complex  potential  w  »  4>  +  iiji,  where  <J> 
is  the  potential  function  and  41  the  Btream  function,  can  be  written  as 

w  -  wp  +  wy  +  wR  (10) 

where  wp  represents  the  parallel  flow  around  the  body,  Wy  is  the  contribution  due 
to  the  presence  of  m  vortices,  and  wR  is  the  term  which  takes  account  of  the 
rotation  of  the  body.  Since  the  complex  potential  io  by  definition  the  same  in  the 
a-  and  C-planes,  the  contribution  due  to  parallel  flow  is13 

w  -  U((;e"lu  +  sl  eiut)  (11) 

P  5 

with  U  the  constant  velocity  of  the  parallel  flow  and  a  the  angle  of  attack 
measured  counterclockwise  from  the  positive  real  axis.  The  term  Wy  representing 
the  m  vorticeB  is  given  by 

m  2 

Wy  -  i  1  k.  [log(t  -  -  log (C  -  — )  1  (12) 

^  k»l  *>k 
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where  teR  denotes  the  strength  of  the  kth  vortex  and  its  position.  is  the 
complex  conjugate  of  tR. 

The  contribution  wR  can  be  obtained  with  a  method  described  by  Milne- 
Thomson.13  For  a  body  rotating  with  angular  velocity  to(t),  the  stream  function 
at  a  point  o  of  the  circle  is  given  by 

2 

2iip  ■  iwss-  iw  f(c)  f(— )  »  iuiB(o)  (13) 

o 

where  B(o)  is  called  the  boundary  function.  If  B(o)  is  written  in  the  form  of  a 
Laurent  series  in  o,  then  B(o)  ■  B^(a)  +  B2(o)  where  B1  contains  the  negative 
powers  of  o  and  B2  the  non-negative  powers.  With  the  aid  of  Cauchy's  residue 
theorem13  it  follows  that 


wR<0  ■  iwB^;) 


(14) 


Bryson12  has  discussed  B^(t)  for  the  cane  covered  by  Equation  (3).  The  general 
situation  a  |l  0  requires  the  evaluation  of  a  double  sum  which  converges  slowly  near 
the  body  and  more  rapidly  farther  away  (sec  Appendix  B) .  For  a  ■  0  this  double 
sum  can  be  reduced  according  to  Bryson  12  to 


I 


H+l 


"Cr 


r(l  +£)enS,+2  -f-- 


r(*+-+i) 


.  2ir 

sin 


(a-0) 


(15) 


where  F  is  the  gamma  function.  This  expression  is  undetermined  for  n 
However,  one  easily  obtains  from  Equation  (3): 
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c 


for  n«l,  a»0 
for  n-2,  a"0 


1  and  2. 


(16) 

(17) 


The  angle  of  attack  u  is  related  to  the  angular  velocity  at  by 

t 

a  ■  an  +  /  w  dt  (IB) 

0 

The  complex  conjugate  velocity  of  the  flow  is  given,  except  for  the  point 
vortices  themselveB,  by: 


dz  .  dw  dw  dt 

— -■U-iV»-T““-TTj 
dt  dz  dt  dz 


(19) 
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where 
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•ia) 
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‘"fc 
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2>+iwd f 


(20) 
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ds 
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-  48 
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-  A" 


-1 


“1 


(21) 


The  velocity  of  the  k  vortex  le,  according  to  Routh'e  theorem  (Milns-Thomson13) , 


with 


A  i 

dt  dt  d«k  “  2  Kk  d;k 


+  1 


A. 

5"ck 


(22) 


When  the  reference  frame  is  to  be  fixed  to  the  body,  the  solid-body  rotation 

w  ■  -  -|iuC  must  be  added  to  the  terms  in  Equation  (10).  The  corresponding 
velocity  is 

ft-iw?  (23) 

3.2  DIMENSIONLESS  FORM  OF  THE  COMPLEX  VELOCITY 

The  complex  velocity  of  the  flow  field,  including  that  of  the  vortices, 
Equations  (19)  and  (22),  is  made  dimensionless  by  U.  For  the  coordinates  and  other 
quantities  with  the  dimension  of  length,  the  parameter  2b  is  chosen  as  the  charac¬ 
teristic  length.  In  particular,  the  time  t  1b  made  dimensionless  by  U/2b  and  the 
vortex  strength  k  by  l/2b  U.  The  Bame  notation  is  used  for  the  nondimensional 
quantities  in  the  resulta  (Section  6). 
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A.  VORTEX  GENERATION  AND  SHEDDING 


A . 1  FEEDING  MECHANISM 

The  process  of  vortex  generation  requires  the  existence  of  a  boundary  layer 

and  its  separation  from  the  body.  After  separation  the  boundary  layer  becomes  a 

free  shear  layer  in  which  vorticity  ie  concentrated  in  a  thin  layer  for  high 

Reynolds  numbers,  Such  a  layer  rolls  up  In  time  by  forming  a  vortex.  In  potential 

flowt  which  does  not  allow  the  creation  and  elimination  of  vorticity,  the  roll-up 

is  modeled  by  a  discontinuity  line.  Self-similarity  of  the  discontinuity  spiral 

with  advancing  time  is  assumed  for  the  initial  period  whan  the  finitenesa  of  the 

plate  is  not  yet  felt.  This  idea  goes  back  to  Prandtl.14  More  recently,  detailed 

Btudies  were  made  by  Wedemeyur , 1 5  Blendermann, 16  and  Pullin.17  The  discontinuity 

*•2/3 

spiral  close  to  the  center  has  the  form  r  %  $  '  and  is,  therefore,  a  hyperbolic 

type  of  spiral  with  Infinitely  many  turns. 

Numerically,  the  discontinuity  line  itBelf  is  often  approximated  by  a  row  of 
point  vortices  (discrete-vortex  model).  Earlier  difficulties  with  such  a  model 
have  been  largely  overcome  by  a  rediscretization  procedure  developed  by  Fink  and 
Soh.1  This  method  is  used  in  the  present  work. 

At  the  sharp  edge  of  a  body  flow  always  separates  when  it  meets  the  body  under 
a  nonzero  angle.  The  discontinuity  line  originates  at  the  sharp  edge,  It  grows  in 
time  with  new  line  elements  forming  at  the  edge.  In  u  discrete-vortex  model  new 
point  vortices  are  introduced  after  the  time  interval  At.  The  feeding  mechanism 
works  In  the  following  way: 

(a)  At  timo  t+At  vortices  are  convected  away  from  the  edge.  Their  new  positions 
ure  computed  with  the  aid  of  Equation  (22). 

(b)  Each  discontinuity  line  is  rediscretizcd  in  the  physical  plane  so  that  every 

vortex  on  it  lies  at  the  center  of  the  segment  represented  by  the  vortex.  If  is 

the  distance  of  vortex  k  along  the  line  from  the  firBt  vortex  on  the  line,  then 

the  total  length  S  of  the  line  is  bu,  where  M  is  the  number  of  vortices  on  the  line. 
The  segment  length  after  rediscretization  is  given  by  As  ■  S/(M-1)  and  so  the  new 
positions  of  the  vortices  can  be  calculated  from 

sk  -  (k-l)As  (2A) 

where  s,  Is  the  distance  arter  rediscretization  of  vortex  k  from  the  first  vortex 
k 

on  the  line.  The  positions  of  the  first  and  last  vortices  are  unchanged  by  the 
rediscretization  procedure. 
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(c)  The  strengths  of  the  vortices  are  recalculated  according  to  Shouff  10  to 
account  for  the  changed  positions  of  the  vortices.  First,  the  strength  density 
near  each  vortex  before  rediscretization  is  computed  as 


k1/<S2  “  "l* 

if 

k-1 

V(,k+i  “  Vi> 

if 

l<k<M 

~  *M-1* 

if 

k-M 

(25) 


For  the  redistributed  vortices,  the  strength  density  is  approximated  by 


Vl  +  Vl*  r 


1-1 


’l- 1 


(26) 


where  1  has  been  determined  so  that  a^  5  <  a^.  Thus,  the  strength!  of  the 

redistributed  vortices  are  y^As  to  a  first  approximation.  But  Bince  the  procedure 
outlined  so  far  does  not  necessarily  conserve  the  total  vortex  strength  in  each 
discontinuity  line,  the  deficit  or  excess  strength  is  removed  by  adding  an  equal 
amount  of  strength  to  each  vortex.  Hence,  the  new  strengths  are  given  by 

M  M 

K.  -  y.  As  +  -  (  £  k,  -  Z  y8As)  (27) 

k  k  M  1-1  *  1-1  * 


(d)  In  each  discontinuity  line  a  new  vortex  is  introduced  between  the  edge  and  the 
first  vortex  at  a  point  1/3  of  the  distance  from  the  edge  in  the  physical  plane. 

The  c-plane  positions  of  all  the  vortices  are  calculated.  Then  the  strengths  of 
the  nascent  vortices  are  determined  by  satisfying  the  Kutta-Joukowaky  condition 


(&) 


-  0 


1 


(28) 


at  each  fin  tip  Cj  in  the  C-plane.  For  n  fins  n  linear  equations  of  Equation  (28) 
type  must  be  Bolved. 


Accuracy  is  checked  by  computing  the  shape  and  position  of  the  vortex  spiral 
and  the  increase  of  the  total  vortex  strength  at  the  i11'  fin  with  time  for 
various  At t 


KL(t) 


M 

Z  K 
k-1 


ik 


(29) 
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Figure  4a  shows  that  with  decreasing  At  the  number  of  the  Inner  loops  increases , 
but  this  does  not  seriously  effect  the  shape  and  location  of  the  spiral  or  its 
strength  (see  Section  6  and  Figures  4a,  b,  c) . 

4.2  INITIAL  CONDITION 

Starting  a  vortex  sheet  at  t«0  in  a  potential  flow  of  constant  velocity  U 
corresponds  to  the  abrupt  start  of  the  body  from  rest  to  the  velocity  -U.  The 
initial  sheet  for  t  ■  At  can  be  taken  from  the  self-similar  solution  for  a  vortex 
spiral  behind  the  edge  of  a  semi-infinite  plate.  A  trial-and-error  approach, 
however,  shows  that  the  development  of  the  spiral  row  of  point  vortices  is  quite 
insensitive  to  the  placement  of  the  first  vortex  with  respect  to  the  subsequent 
roll-up.  Tha  strength  of  the  first  vortex  again  is  determined  by  the  Kutta- 
Joukowsky  condition  (28).  The  location  of  the  first  vortex,  which  is  arbitrary, 
can  be  described  by  the  distance  As  away  from  the  tip  and  tha  angle  8  between  the 
extension  of  the  fin  and  tho  line  drawn  from  the  tip  to  the  vortex.  Although 
variation  in  As  (from  As  *  '•  J25  to  0.03)  was  not  noticeable  in  the  results,  small 
but  still  insignificant  differences  occurred  when  3  was  varied  (Figure  3). 

4.3  CUT-OFF  PROCEDURES 

The  infinite  turns  in  the  vortex  splrul  cunnot  be  represented  by  a  row  of  point 
vortices,  and  these  turns  are  physically  unrealistic  anyway  (Section  4.4),  Some¬ 
where  the  spiral  has  to  be  cut  off.  Investigations  by  Wedemeyer’&  and  Pullin’ ; 
have  revealed  that  the  almost  circular  windings,  which  represent  the  core  of  the 
vortex,  can  be  replaced  by  u  single  vortex.  Even  this  single  vortex  uppeurs  not  to 
be  necessary  (Fink  and  Soh1).  The  overall  solution  Is  quite  insensitive  to 
arbitrary  cut-off. 

Another  cut-off  procedure  Is  necessury  when  the  rolled-up  spiral  separates  from 
the  body  and  becomes  a  detached  vortex  which  swims  away  in  the  wake.  Although  the 
development  of  a  vortex  row  without  the  use  of  'he  rediscretization  technique  some¬ 
how  takes  cure  of  this  separation  by  itself  (hcc  Figure  7b,  page  21),  the  line  of 
redistributed  vortices  has  to  be  severed  by  a  proper  criterion.  Shoaff11’  uses  the 
condition  dic/dt  ■  min  after  a  certain  developing  time  of  the  vortex  row.  This 
criterion  has  been  applied  here  with  varying  success  (see  Section  6).  The 
detached  line  is,  by  the  way,  also  redlscretixed ,  Shoaff's  technique  of  replacing 
the  detached  line  after  two  body  lengths  by  a  single  vortex  has  also  been  adopted. 


Figure  4  -  Accuracy  Tests  with  Flows  Fast  a  Flat  Plate  at 

a  •  90° 


- -  0.10 

• - At  ■  0.00 

— - At  *  0.020 

- At  -  0.01 


Figure  4a  -  Discontinuity  Lines  at  t  ■  0,5  for  Various  At 
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Figure  4  {Continued) 


Sequence  of  Discontinuity  Lines  at  Selected  Times  for 
Various  fit 


t-  2.0 


- 0 


0  MB 
1  46  d«B 
90  Mg 


Figure  5  -  Accuracy  Teats  with  Flows  Past  a  Flat  Plate  at 
a  ■  90°  for  Various  Initial  Conditions.  Discontinuity 
Lines  at  t  ■  2  with  At  ■  0.025  for  Three  Different  p's. 


.‘.wauuLu.’juii  aiAfuLvii;  a: 
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4.4  VISCOUS  EFFECTS 

The  neglect  of  viscosity  may  restrict  the  usefulness  of  the  model.  There  are 
three  major  flow  regions  in  which  viscosity  cannot  be  ignored,  and  the  influence 
of  such  neglect  on  the  overall  flow  characteristics  must  be  questioned. 

Infinitely  many  turns  of  the  vortex  spiral  are  obviously  not  realistic,  since 
during  the  creation  of  the  vortex  the  core  is  in  a  state  of  solid-body  rotation. 
Moore  and  Saffman19  have  discussed  the  structure  of  the  vortex  spiral  and  have 
given  estimates  of  the  viscous  core  which  is  present  from  t*0  on.  After  separation 
from  the  body  the  vortex  decays  through  dissipation,  and  this  effect  is  not 
simulated  in  the  point-vortex  model  either. 

Vortices  or  blobs  of  vorticity  of  opposite  sign,  which  approach  each  other,  are 
eliminated  or  coalesce  in  a  viscous  fluid.  Also,  when  vortices  approach  a  solid 
surface,  they  are  weakened  or  destroyed  by  the  opposite  vorticity  produced  at  the 
wall.  None  of  these  effects  are  simulated  in  the  inviscid-f low  model. 

4.3  OTHER  SEPARATION  POINTS 

So  far,  flow  separation  has  been  considered  only  at  the  tips  of  the  fins. 
However,  other  separation  points  at  the  body  surface  may  occur  as,  for  instance, 
in  the  case  of  the  circular  cylinder  with  one  fin.  Here,  a  separation  point  must 
exist  on  the  side  of  the  body  opposite  the  single  fin.  Separated  regions  may  also 
appear  between  fins  when  n>  l, 

The  occurrence  of  separation  points  can  be  predicted  with  the  aid  of  boundary- 
layer  theory.  Such  a  prediction  method,  together  with  a  technique  to  provide  for  a 
vortex  sheet  at  the  separation  point,  will  be  included  in  the  computer  code  at  a 
later  time. 


5.  FLOW  CHART 

The  sequence  of  computations  ia  indicated  in  the  following  flow  churt.  The 
calculation  of  force  coefficients  is  included. 

6.  SOME  RESULTS 

A  computer  program  based  on  the  equations  of  Sections  2  through  S  has  been 
developed  for  circular  cylinders  under  the  restriction  n  £  4.  However,  most  cases 
of  practical  interest  are  covered  under  this  restriction. 
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COMPUTE  VELOCITIES  OF  THE  VORTICES 
IN  COMPUTATIONAL  AND  PHYSICAL  PLANES 
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Two  samples  have  been  selected  for  comparing  resultB  of  the  present  computer 

program  with  those  from  the  literature.  A  third  example  gives  new  results.  More 
complicated  cases  of  practical  interest  will  be  published  later. 

6.1  VORTEX  SHEDDING  FROM  A  FLAT  PLATE  AT  a  -  90* 

This  case  is  of  particular  interest  aince  its  results  can  be  compared  with 
results  of  Fink  and  Soh,1  which  are  based  on  a  similar  point-vortax  method,  and 
with  analytical  results  of  Wedemeyer15  for  a  discontinuity  line,  at  least  for  the 
initial  phase  of  the  toll-up.  The  results  are  presented  in  Figures  Ac  and  6. 

In  Figure  6  the  roll-up  of  the  discontinuity  line  is  shown  at  times  t  -  0.25, 
0.5,  0.75,  1.0,  and  1.5  with  At  *»  0.02.  Up  to  t  ■  0.75  the  curves  are  compared 
with  Wedomeyer's  self-similar  solution,15  which  is  valid  for  a  semi-infinite  plate. 
From  t  ■  0.75  on,  deviations  occur  because  of  the  influence  of  the  finltaness  oi 
the  plate  width.  (According  to  Wedemeyer15  differences  between  lnfiniteLy  wide 
and  finite-width  plates  become  noticeable  from  t  ■  0.6  on.)  For  t  *  1.0  and  1.5 
the  results  are  compared  with  those  of  Fink  and  Soh.1  The  agreement  in  quite  good. 
In  all  cases  the  limiting  curve  for  t  ■  »  by  Helmholtz  is  shown,  along  which  the 
discontinuity  Bpiral  rolls  up  until  it  becomes  unstable.  It  may  be  mentioned  that 
in  this  early  phase  the  flow  is  symmetric,  and  no  attempt  has  been  made  to  induce 
alternating  vortex  shedding  through  an  initial  asymmetric  disturbance. 

In  Figure  Ac  the  Increase  in  total  vortex  strength  k  with  time  for  one  half  of 
the  plate  is  compared  with  the  corresponding  result  by  Fink  and  Soh.1  Their  data 
are  slightly  larger  and  agree  with  those  of  Wedemeyer  for  a  plate  of  finite  width. 

6.2  VORTEX  SHEDDING  FROM  A  FLAT  PLATE  AT  a  -  A5° 

Results  for  this  case  can  also  be  compared  with  those  in  the  literature. 

In  Figure  7  the  roll-up  of  the  discontinuity  line  is  shown  for  t  *  0.5,  1.0,  and 
2.0  with  At  »  0.05  and  is  compared  with  the  curves  by  Belotserkovskii  and  Nisht,20 
who  did  not  use  a  redisc.retization  procedure.  The  advantage  of  rediscretization  is 
particularly  demonstrated  for  the  roll-up  behind  the  leading  edge.  In  Figure  8 
the  same  situation  is  presented  at  the  slightly  different  times  t  *  0.39,  1.12,  and 
1.87  and  the  results  are  compared  with  the  corresponding  solutions  of  the  Navier- 
Stokes  equations  for  Re  ■  2bU/v  ■  200.  The  flat  plate  1b  here  approximated  by  a 
thin  elliptic  cylinder  of  infinite  length  with  a  width-to-thickness  ratio  of  10  to 
1.  A  discussion  of  the  differences  has  already  been  published.21 
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Figure  7  -  Development  of  Discontinuity  Lines  Behind  a  Flat  Plate 


In  Figure  9a  the  total  strength  of  the  leading-edge  vortex  is  compared 
with  that  of  the  trailing-edge  vortex  .  The  absolute  amount  of  the  letter  ia 
■lightly  smaller  so  that  the  sum  of  the  two  is  not  aero  (Figure  9b) .  This  violates 
the  conservation  law  of  vorticity.  In  reality ,  however,  boundary  layers  contribute 
to  the  generation  of  vorticity  which  would  account  for  the  difference.  In  an 
inviacid-flow  model  a  bound  vortex  cau  be  introduced  at  the  center  of  the  circle 
to  balance  the  difference.  This  technique  was  studied  with  the  present  computer 
program.  Since  the  results  were  not  significantly  different,  the  incorporation  of 
auch  a  bound  vortex  was  abandoned. 

Difficulties  have  been  encountered  with  the  cut-off  procedure.  The  trailing- 
edge  vortex  could  be  separated  with  Shoaff's  criterion  (Section  4.3)  and  could 
simulate  vortex  shedding  satisfactorily.  However,  the  leeding-edge  vortex  did  not 
move  away  fast  enough  after  cut-off  and  intarfered  with  the  development  of  the  new 
vortex  sheet.  The  problem  has  not  yet  been  solved,  but  it  may  be  a  consequence  of 
the  rediscretization  procedure. 

6.3  VORTEX  SHEDDING  FROM  A  CIRCULAR  CYLINDER  WITH  TWO 

FINS  AT  a  •  45* 

In  the  final  example  the  vortex  shedding  from  a  circular  cylinder  with  two  fins 
at  a  -  45#  was  compared  with  that  from  a  fiat  plate  at  a  ■  45*.  Figure  10  displays 
the  development  of  the  discontinuity  lines  for  both  cases.  Up  to  t  ■  0.6  the 
leading-edge  vortices  do  not  show  any  visible  differences,  but  the  trailing-edge 
vortex  for  the  flat  plate  ia  slightly  stronger.  Beyond  t  -  0.6  the  leading-edRe 
vortex  is  deformed  by  the  presence  uf  the  cylinder.  The  corresponding  data  for  the 
total  strengths  and  K,,  in  Figure  9  confirm  that  slightly  smaller  than  K,, 

for  the  flat  plate. 

The  two-dimensional  time  development  of  the  discontinuity  line  in  Figure  10  can 
also  be  interpreted  bb  a  special  growth  in  a  three-dimensional  flow  within  the 
frame-work  of  a  strip  theory.  Then,  t  iB  replaced  by  the  coordinate  Z  perpendiculur 
to  the  x,y-plane  with  the  aid  of  the  constant  velocity  W  in  the  Z-directions 
t  ■  Z/W.  A  computer-generated  perspective  view  is  presented  in  Figure  11. 


Figure  9  -  Total  Strengths  of  the  Vortices  at  the  Trading  and  Trailing 
Edges  for  Flows  Past  a  Flat  Plate  and  Past  a  Circular  Cylinder 
with  Two  Fins  at  a  =  45® 
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7.  CONCLUSIONS  AND  SUMMARY 


a.  A  computer  program  haa  bean  developed  to  simulate,  by  means  of  a  discrete- 
vortex  model,  vortex  shedding  from  a  circular  cylinder  with  up  to  four  evenly 
distributed  fins. 

b.  A  new  approximate  conformal  transformation  for  a  circular  cylinder  with  two 
dihedral  fins  (Figure  2)  has  been  derived,  and  a  conformal  transformation  for  an 
elliptic  cylinder  with  two  fins  haa  been  improved  (Section  2). 

c.  A  numerical  method  haa  been  devised  for  computing  the  coefficients  of  the 
Laurent  aeries,  which  are  necessary  to  find  the  potential  function  for  the  rota¬ 
tional  motion  of  the  body  (Section  3.1  and  Appendix  B) . 

d.  The  feeding  mechanism,  that  is,  the  Introduction  of  the  new  point  vortex  near 

the  tip  at  each  time  step,  is  the  crucial  process  in  the  whole  model.  Although  the 

shape  and  the  location  of  the  spiral  discontinuity  lines  are  quite  Insensitive  to 
various  approximation  schemes,  the  strengths  of  the  discrete  vortex  rows  are  some¬ 
what  weaker  than  those  reported  by  Wedemeyer*5  and  Fink  and  Soh.1 

e.  Although  results  for  the  force  coefficients  are  not  reported  hers,  preliminary 
studies  indicate  that  they  are  very  sensitive  to  the  kind  of  feeding  mechanism 
used.  This  is  also  reflected  in  the  different  results  of  Belotserkovskii  and 
Nisht20  and  Barpkaya22  for  the  vortex  shedding  from  an  inclined  plate. 

f.  The  CP  time  in  seconds  cm  the  TI-ASC  is  equal  to  0.003  for  bodies  with  two 
fins. 

8.  PROPOSED  EXTENSIONS  AND  REFINEMENTS  TO  THE  PROGRAM 

The  usefulness  of  the  computer  program  described  can  be  enhanced  by  Incorpora¬ 

ting  the  following  extensions  and  refinements: 

a.  Include  the  roll-up  of  discontinuity  lines  shed  from  the  cylinder  (other  than 
the  tips  of  the  fins).  TUIb  requires  building  in  a  boundary-layer  code  for  de¬ 
termining  the  point  of  separation  and  the  amount  of  vorticity  shed  at  that  point. 

b.  Improve  the  feeding  mechanism  to  obtuin  reliable  force  and  moment  coefficients. 

c.  Include  the  computation  of  force  and  moment  coefficients. 

d.  Investigate  viscous  effects  and  make  appropriate  corrections. 

e.  Consider  other  conformal  transformations  of  practical  interest. 
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APPENDIX  A 

DERIVATION  OF  THE  CONFORMAL  MAPPING, 
EQUATIONS  (4)  AND  (5) 


The  derivation  of  the  conformal  mapping  of  the  finned  cylinder  in  Figure  2  onto 
a  circle  is  divided  into  the  following  steps t 

(1)  Conformal  transformation  of  the  original  figure  in  the  z-plane  onto  the 
auxiliary  plane  z^,  in  which  the  fins  are  mapped  to  arcs  of  a  hyperbola  (Figure  12): 


z.  -  z  +  — 
1  z 


Then,  the  points  A  through  H  and  the  corresponding  points  A^  through  in 
Figure  12  are  given  by 


A  ■  +  a 


B  -  -  ae 


D  -  B 


Ai  -  +  2a 

IL  ■  **  2a  cos  9 
-2 


C^  ■  -  (b  +  ^-)cos  9  +  l(b  -  ^-)ain  6 


D1  "  B1 


F  -  -  aeJ 


H  -  F 


F1"B1 

2  2 
-  -  (b  +  ^-)cos  0  -  i(b  -  ^-)sin  0 

Hl"  B1 


The  lengths  k  and  l  are 

k  -  -•  (CL  -  C1)  «  |  (b  -  *~)  sin  0  (A. 

II  -  Dj_  -  |  (C  +G1)  «  ^  (b-a)2  cos  0  (A1 

The  hyperbolic  arc  1“  now  approximated  by  the  circular  arc  through  these 

points  with  the  radius 

R  -  ~  (4k2  +  t2)  ^ A! 

(2)  According  to  Betz23  this  circular  arc  can  be  mapped  onto  a  circle  by  means  of 


,  -  —  -  z.  +  (b  +  — )  cos  9 

2  Zj  1  b 
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Figure  12  -  Sequence  of  Conformal  Mappings  from  a  Circle  with  Two 
Dihedral  Fins  onto  a  Circle 


with 


.2  2 

A0  "  positive  root  of  z,  -  -  2a  +  (b  +  4-)  cos  8 

■  •  2 «  b 


B2  ■  |  “  /4k2  +  t2  -  |  +  f\ B  ~l 


(A7> 


•tc., 


where 


M  -  ,  r  -  i(B,  -  D-)  -  /H7 


2'  2  2  2 


(A8) 


(3)  Now  the  figure  in  the  mapped  onto  a  straight  line  in  the  z^-plane 

by  2 

z3  “  z2  "  2  +  z 2  -  i/2  (A9) 

All  points  of  the  figure  lie  on  the  real  axis  with 

A3  “  A2  “  2  +  A^-T/2 


Ba  “  Ba  “  1  +  ifriTi  -  ^1? 


(A10) 


etc. 

(4)  Finally,  the  straight  line  in  the  z^-plane  Is  mapped  onto  the  circle  in  the 
C-plane  with  radiuB  c  by 

2 

5  +  ■  z3  -  »  (All) 

with 

8  -  ~  (A3  +  E3)  (A12) 

c  -  \  (A3  -  E3)  (A13) 

(See  Figure  12.)  Combining  the  four  transformations  to  two,  one  arrives  at 
Equations  (4)  and  (5). 
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APPENDIX  B 

LAURENT  SERIES  FOR  w„ 


The  conformal  mapping  z  ■  f(C)  is  expreseed  in  an  infinite  series  of  the  form 


c  3 

z  ■  I  a  (r)  »  where  a  ,  -  c 

q— 1  q  4 


(Bl) 


The  principal  part  of  zz  on  the  circle  a  (see  Equation  (13)),  which  is  required  to 
obtain  B1(?) ,  is 


P.P.(zz)  »  E  b  (~)q  on  Ul  “  n 
q-1  q 


(B2) 


Since  It  is  very  laborious  to  determine  a.  and  b  analytically,  these  coefficients 

q  q 

are  computed  numerically.  From  Equation  (Bl)  it  follows  that  on  the  circle 


z  ■  E  a  e 
q*-l  q 


-iq9 


ie 


where  t  ■  ce  .  Equation  (B3)  is  a  Fourier  series  whose  coefficients  can  be 
determined  from 

2  it 


a  I  « •l<'9  de 

1  2"  0 


or  in  discretized  form  from 


AO 


N 


iqe, 


11  ^  nfr  E  2.  e 

q  2%  e-1  i 


The  coefficients  a  are  real  if 

q 


because 


z(0)  -  z(2ti-9) 


a  -■—•[/*  elq3dO  +  /  i  eiqU  dB] 
q  2*  0  n 


-  ~  J  (z  elq0  +  z  s’lq9)dG 
■  —  /  Re(ze.^q0)dO 


(B3) 


<B4) 


(B5) 


(Hfi) 


(B7) 
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Similarly,  Equation  (B2)  becomes  on  the  circle 


■i 

,1 


whose  coefficients  b  are 


P.P.(«i)  -  l  be'1’9 
<1-1  ’ 


\  ■  n  •1”  de 


*  bl  "  ,  -  .‘""i 


<B8) 


(B9) 


Again,  for  the  symmetry  condition  (B6)  It  can  be  shown  In  a  way  analogous  to 
Equation  (B7)  that  b^  is  real. 

As  an  example,  for  a  circular  cylinder  with  two  fins  (n-2)  and  varioua  ratios 
a/b  the  coefficients  a^  and  b^  are  displayed  in  Tables  1  and  2  for  q«l  through  50 
with  A6  •  2ir/1080.  For  a/b  “0,  a^,  ■  b_  ■  0  except  a,  ■  0.25,  b.  ■  c^  ■  0.0625. 

C]  ()  X  X 

For  the  other  extremum  a/b  ■  1,  a  *  ;  and  thuo  all  coefficients  are  aero.  Despite 
the  fineness  of  A0  the  data  in  both  tables  are  accurate  only  to  i  5  -  06  because  of 
round-off  errors. 
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TABLE  1  -  COEFFICIENTS  a  FOR  n  -  2, 

q 

a:b  ■  0.2,  0.4,  0.6,  0.8 


n 

wmm 

0.4 

0.6 

O.B 

i 

.22152*00 

.15216*00 

. 75256-01 

.15766-01 

3 

.32776-01 

.72326-01 

.56616-01 

.16*06-01 

5 

-.23071-01 

-.35156-02 

,32666-ni 

.16556- Cl 

7 

.12112-01 

-.17666-01 

,80516-02 

.14456-01 

0 

-.27112-02 

.2*236-"? 

-.6*4  *6-P2 

•  11596-01 

11 

-.31526-02 

.6735E-12 

-.1M4F-32 

.6**56-02 

13 

.50806-02 

-.2165F-02 

• ,462  4E -0  2 

,6-*i6r-02 

15 

-.38166-02 

-  ,62*66-"2 

» 1*  6  4F  -0  2 

.25706-02 

1? 

.10516-02 

.18716-32 

.461 36-02 

»02*f6- 03 

15 

.1 3646-02 

,3*476-02 

.*2*66  -02 

-.12556-02 

21 

-.24526-02 

-.16576-02 

-.7664E-0? 

-. 22636-02 

23 

.15576-02 

-.25156-02 

-.375  *F «P  2 

-.26046-02 

2? 

-.60216-03 

.14O0F-C2 

-.247  76-0  2 

-.24226- 02 

2? 

-.81156-03 

» 1  64*6- C  2 

-.?C>Mr-P3 

-.1 *6?E- 02 

2? 

.15106-1*2 

-.1**36- :? 

.171  ?f-o  2 

-,10  526-  02 

31 

-.1267F-02 

-  .  1 36  ?*• *  2 

,15476-r  ? 

-.27606-03 

33 

.3*306-03 

. 12366-  0  2 

,66  Ptf «c  j 

.44426-  03 

e m 

,56if F-P3 

,10426-'? 

-.10416-02 

,4*556- 03 

KB 

-.1(506-02 

-.  113  36-  "2 

-.15476-02 

.1  2336-  02 

im 

,««CCF-"3 

-.76406-"! 

-.*1056-0? 

.12406-02 

41 

-.26416-03 

.  10426-32 

,  R5  C  46  •  0  3 

.1*256-  I)? 

63 

-.421*6-05 

,•3266-?? 

.1*266-02 

.66466-03 

45 

.**456-0! 

-  *  «m  «r-  o  3 

,7bu*r-o? 

.21626-0? 

47 

-.66636-03 

-,42?8F-"3 

».?73*F«0* 

-.207  66-03 

40 

.10166-?! 

•  *»i  r.f- "  * 

-.qq^uF-o? 

-»64?2r- 03 

2 


.oJf-t-G  1 

•  Ja24e-02 
«  32  33&  -g  2 
«22.1e  &  -  0  2 

•  lit  7t-02 

•  loUofc  -03 

2  b  J  c.  —  J  3 

•  »459e-Ci  3 
« -3-0e  -03 
.12426.-03 

•  le»0-fi-04 

•  2(»7bt  -04 
.23976-03 
.61226-0- 
.4527  t-U- 
.1/  -*.6-03 
.15556-03 
.pi  -7e  -0- 

•  5- 0  It -0- 
.1H.U- 03 
. 1 1 uec -03 
.  -3-J 1  -0"* 
•36026-0- 
.  t»7  9*  t  -0  - 
. rt35lt -«<•* 


0.8 


5-4  7£-ol 
l*47r.-0  1 
3»3oh-0  2 
34«j3c-0  2 
ho;e-o2 

ln^ot-u  2 

12 936-0 2 

44446-03 
•92756-u  3 
5127r.-03 
.72006-03 

•  30-7<i-C3 
.  i79-t-C3 
.177  OI--C3 
.•.7792-03 

•  96  0«‘t-C4 
3*-#  <6-g  3 

.  3*  i  Sr-04 
.  33  1« -0  3 
.  4i  23t-l  to 


32261-01 
22  9-6-01 
10  6  -It-  0  1 
1- 3  66*  02 
32  w lt-02 
32  7 6t- 02 
*6416-03 
11716-02 
170-6-02 
77066-03 
00 u76- 03 
10556-02 
0-436- 03 
14106-03 
4  0  256-03 
»-636* 0J 
CS 7*6-0- 
-7136-03 
-0656- 03 
73  5  ft- C- 


32-16-03 
39  .i  C6-  w  3 
131-6-03 
21 1-6-03 
33516-03 


. 93956-02 
. 87-2  E-C  2 
.  77216-02 
*  6-2*E-02 
. -9746-02 
.3-956-02 
.  20  99  E  -  02 
. 64696-03 
•»  bio5t-04 
-.71426-03 
-.  10o<SE-02 
-. 11456-02 
-.  10u7t-02 
-. 72296-03 
-. 35906-03 
-.  Io926-04 
.  27  53  6-  03 
. ho* 16  -  0  3 
. 5-776-03 
. 51576-03 
.39626-03 
. 21 *76  -  03 
.  27  9.36  -0- 
*• 14536-03 
•»  27ur»t-  03 
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DTNSRDC  ISSUES  THREE  TYPES  0?  REPORTS 

.Pi,  DTNSRDC  REPORTS,. A  FORMAL  SERIES,  CONTAIN  INFORMATION  Of  PERMANENT  'r 
NWAU  VALUE.  THEV  WRRY^A  CONSECUTIVE  NUMERICAL  !0£NTIFICA7 SON  REGARDLESS 
7HEUI  0USS!R!CAT!ONO'R  THE  ORIGINATING  D*>AR»T*tNT. 

DEPARTMENTAL  REPORTS,  A  ASMS  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  «»f 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  Of  U’AJTEt?  INTEREST  OP  SlGNiFiC.V 
THEY  CARRY  a  DEPARTMENTAL  Ai.fTtANUMEP.iOAL  ICRNTlFICATiON. 

•3.  TECHNICAL  MEMORANDA  AN.  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DC  Cm  1  h  >  A 1 
Q?  LIMITED  USE  AND  INTEREST,  THEY  ARE  PRIMARILY  WORKING.  PAPERS  IW  LNDSO  FT - 
TERN’AL  USE.  THEY  CARRY  AN  IDENTIFYING  fTJftfSSS  WHICH  INDICATES  THStfi  TYPE  ANP 
•NUMERICAL  CODE  of -'THE  ORIulNATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  D  r.\  " 
MUST  8E  APPROVED  BY  THE  HEAD  QF  TnE  ORIGINATING  DEPARTMENT  ON  A  CASS-O.v-0 
BASIS. 


